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Abstract:  In  this  paper,  singularly  perturbed  nonlinear  differential/algebraic  equations 
(DAEs)  are  considered  and  a  proof  of  the  existence  and  uniqueness  of  a  solution  is  given. 
Asymptotic  expansions  for  such  a  solution  are  obtained  and  proved  to  be  uniformly  conver¬ 
gent.  This  generalizes  known  results  about  asymptotic  expansions  of  singularly  perturbed 
ordinary  differential  equations. 

AMS  subject  classification:  41A60,  34D15,  34E05 


1.  Introduction 

We  consider  the  asymptotic  behavior  of  solutions  of  a  singularly  perturbed  DAE  of  the 

form 

=  /i(z,y,2.e) 

«/'  =  f2(x,y,z,e)  (1-la) 

0  =  f3(x.,y,z,e) 

together  with  the  initial  conditions 

*(0,e)  =  f(e),  y( 0,e)  =  ij(c),  *(0,e)  =  <(e)  (1.1b) 

where  x  €  Rm,  y  €  Rn,  z  G  Rk  and  f.  €  jR1.  It  might  be  expected  that  the  theorems 
developed  for  singularly  perturbed  ODEs  can  be  used  to  study  the  singularly  perturbed  DAE 
(1.1).  For  this  it  would  be  natural  to  apply  a  standard  index  reduction  to  (l.la,b);  that  is 
to  reduce  (l.la,b)  to  a  singularly  perturbed  ODE  by  differentiating  the  constraint  equation 
f3(x(t,e),  y(t,  e),  z(t,e),  e)  =  0  with  respect  to  t  along  any  solution  (r(t,e),  y{t,e),  z(t,  e)) 
of  (l.la,b).  Then  by  applying  the  known  theory  for  singularly  perturbed  ODE’s  (see  [Ho] 
or  [OM2]),  we  might  expect  to  obtain  the  desired  results  for  the  singularly  perturbed  DAE 

*  This  work  was  supported  in  part  by  ONR-grant  N-00014-90-J-1025,  and  NSF-grant  CCR- 
9203488. 


1 


(l.la,b).  But  this  is  not  the  case.  In  fact,  by  differentiating  the  constraint  equation  in  thi 
system  (1.1a, b)  with  respect  to  t  along  the  solution  (x(t ,  e),  y(t,  e),  z(f,  e)),  we  obtain 


fy'  =  f2(x,y,z,t), 

ez'  ~  -(Dtf3(x,y,z,e))-l(eDxf3(x,y,z,e)fi(x,y,z,e),  (1.2) 

+  Dyf3{x,y,z,e)f2{x,y,z,e)), 

^(0,  e)  =  $(e),  y(0,e)  =  i|(e),  2(0,  e)  =  <(«). 

For  c  =  0,  the  reduced  problem  for  (1.2)  then  has  the  form: 

*  =  /i(x,y,2,0), 

0  =  /2(x,y,  2,0),  (1.3) 

*(0)  =  f(0). 

This  reduced  system  has  lost  k  constraint  conditions  which  means  that  the  assumptions  of  the 
theorems  in  [Ho]  or  [OM2]  for  singularly  perturbed  ODEs  are  not  satisfied.  Therefore,  we  have 
to  study  the  singularly  perturbed  DAE  (1.1)  directly. 

Under  certain  assumptions,  we  prove  that  (l.lajb)  has  a  unique  solution  on  the  interval 
[0,TJ  for  all  small  e,  for  which  asymptotic  expansions  have  been  obtained  and  proved  to  be 
convergent  uniformly  in  [0,T], 

As  background  for  the  presentation,  Chapter  2  below  presents  a  summary  of  some  known 
existence  results  for  DAE’s,  which  can  be  applied  to  (1.1)  and  its  reduced  system. 

Chapter  3  addresses  the  limiting  problems  in  which  the  reduced  problem  and  the  inner 
and  outer  problems  for  (1.1)  lire  introduced,  and  the  regular  degeration  of  (1.1)  is  defined  as 
well. 

While  asymptotic  expansions  and  existence  of  the  outer  solutions  are  considered  in  Chapter 
4  and  5,  respectively,  inner  solutions  of  (1.1)  are  studied  in  Chapter  6  and  7. 


2.  Background  on  DAE’s 

To  ensure  the  existence  of  solutions  of  (1.1),  wo  impose  the  following  assumption  on  the 
system 

Assumption  (I):  There  are  non-empty  open  sets  Vz  C  Rm ,  T>y  C  Rn,  Vt  C  Rl  and  Jt  C 
R*  i  —  {  e  |  |c|  <  c\  c'  >  0),  such  that  the  mappings  fi:VxJt—>  R"‘ ,  /2  :  V  x  J(  — * *  Rn, 
fz  :  T>  x  J,  Rk  ,  £  :  J(  — »  VT,  t]  :  J,  — ►  Vv  and  (,*  :  J(  —>  T>z,  where  V  =  T>x  x  T>y  x  V{ ,  are 
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continuous  on  the  indicated  domains.  Moreover  for  fixed  e  €  Jt,  /i,  f2,  fz  are  of  class  C°°  on 
V  and  the  initial  point  (£(e),  17(e),  C(e))  satisfies  the  compatibility  condition 


f3(t((),r/(e),<;(e),e)  =  0. 

Finally,  assume  that  f2(z:,y,z,0)  ^  0  and  f3(x,y,z,0)  ^  0  in  the  domain  V,  and  that  the 
Jacobian  matrix 

X>,/s(e o,f?o,<o,0)  (2.1) 

is  nonsingular,  where  Co  =  C(0),  t)q  —  i/(0),  Co  =  C(0)- 


The  infinite  differentiability  of  /],  /2,  fz  is  assumed  here  only  for  the  sake  of  simplicity. 
We  are  interested  in  the  existence  of  solutions  of  (l.la,b)  on  some  interval  [0,T]  where  T  is 
independent  of  e,  and  with  the  asymptotic  behavior  of  the  solutions  of  (l.la,b)  as  t  tends  to 
zero.  For  these  asymptotic  considerations  some  further  conditions  will  be  needed  which  will 
be  stated  in  the  next  subsection. 

For  the  existence  of  solutions  of  (l.la,b)  the  condition  (2.1)  in  Assumption  (I)  ensures  the 
solvability  of  (l.la,b).  Indeed,  from  (2.1)  it  follows  that  there  exists  an  e'  >  0  (d  6  Jt)  such 
that 

is  nonsingular  for  any  fixed  e,  0  <  e  <  e'  .  This  implies  that  the  system  (l.la,b)  is  a  DAE  of 
index  one  in  some  neighborhood  of  (£(e),  »}(e),C(e))- 

The  existence  and  uniqueness  of  a  solution  of  (l.la,b)  will  be  based  on  the  following 
existence  theorem  for  the  solutions  of  initial  value  problem  of  the  form: 

u'  =  Fi(u,v), 

0  =  F2(u,v),  (2.2) 

u(0)  =  uo,  v(0)  =  vo, 


Proposition  2.1:  Suppose  that  the  mappings  F]  :  T>u  x  Vv  C  J?r+a  — ►  Rr  and  F2  :  T>u  x  Vv  C 
Rr+3  — >  R1  are  of  class  C1  on  their  domains  where  Vu  C  Rr  and  T>v  C  Ra  are  non-empty  open 
sets,  and  that  the  initial  point  (uq,Uo)  G  T>u  X  T>v  satisfies 

^(^o^o)  =  0, 

and  the  Jacobian  matrix 


is  nonsingular.  Then  there  exists  a  C1 -solution  (u(t),  t’(f))  £  T>u  x  Vv  of  (2.2)  which  is  unique 
on  some  interval  [0,To],  T0  >  0.  Moreover ,  the  component  u(t )  is  of  class  C 2  on  [0,  To]. 

Proof:  The  result  is  a  consequence  of  the  existence  theorem  for  DAEs  of  the  form  (2.2)  given 
in  see  [Rh4].  i 


By  applying  Proposition  1.1  to  (l.la,b),  we  obtain  the  following  existence  theorem: 

Proposition  2.2:  Under  assumption  (I),  for  any  fixed  e  >  0  (e  <  e'  £  J(),  there  exists  a 
unique  solution  (x(t,e),y(t,e),z(t,  e))  for  the  DAE  ( l.la,h)  on  some  interval  [0,T(],  where  T( 
depends  on  f. 


3.  The  boundary  problems 

In  order  to  study  the  asymptotic  behavior  of  solutions  of  (1.1),  we  formally  set  e  =  0  in 
(1.1a)  and  remove  the  initial  conditions  for  y  and  z,  and  then  obtain  the  system 

Xo  =  MXo,Y0,Z0,0), 

0  =  /2(Xo,r0lZo,0),  (3.1a) 

0  =  /3(.Y0,r0,70,0), 

with  which  we  associate  initial  conditions  of  the  form 

*o(0)  =  «o,  Vo(0)  =  ro°5  Zo(0)  =  Z°.  (3  .It) 

Thus  we  obtain  from  Proposition  1.2  the  existence  result: 

Proposition  3.1:  Under  assumption  (I)  let  {(q,Y$,Z%)  £  V,  £0  =  £(0),  be  a  point  such  that 

M^Y0°,Z°0, 0)  =  0, 

MZoX,Zl0)  =  0 

and  that  the  matrices 

D*A((oX,Z2,0)  (3.2  a) 

B°  =D!/f2((o,Yo0,ZS,()) 

,o  tyL\ 

-  D'MtoXXMiD'moX^ZlO))-'  D9f3(toX>Z°,0) 

are  nonsinguiar.  Then  the  system  (3.1)  has  a  unique  solution  (,Y0(f),  l'o(t),  Z0(t))  €  V  on 
some  interval  [0,T],  which  satisfies  the  initial  condition  (3.1b).  Moreover,  for  t  £  [0,T]  the 
matrices 

DJ3(Xo(t),Yo(t),Zo(t),0),  (3.3a) 
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(3.36) 


B(t)  =D9f2{t0,  K0(<),  Z0(t ),  0)  -  DzMto,Yo(t),  Z0(t),  0) 

(^/3(6,Ko(i),^o(*),0)r1JDy/3(6I3/o(0,2o(<),0) 

are  nonsingular. 

Proof:  Since  the  matrices  (3.2a/b)  are  nonsingular,  there  exists  a  neighborhood  Og  C  V  of 
the  point  such  that  for  all  (x,  y,z)  €  Og  the  matrices 

Dzh(x,y,z,0), 

Dyh(x,y,  z,  0)  -  Dzf2(x,  y,  z ,  0)  (Dzf3(x,  y,  z,  0))-1  Dyf3(x,  y,  z,  0) 

are  nonsingular.  This  implies  that  the  Jacobian  matrix 

(  Dyf2(x,  y,  z,  0)  DJ2(x,y,z,0)\ 

\Dyf3(x,y,z,0 )  Dzf3(x,y,z,0)J 

is  nonsingular  in  Og ,  as  follows  directly  from  the  identity 

(Dyh  DJ2\ 

\Dyh  D.h) 


=  (*m> cm  DzMDzhr^fDyh-DzhiDJ^Dyh  0  N 
V  0  hxk  )\  Dyh  DJ3J- 

Then  a  direct  application  of  the  Proposition  2.1  to  the  system  (3.1a/b)  shows  that  this  system 
has  a  unique  solution  (Xo(t),Y0(t),  Z0{t))  €  Og  e  V  on  some  interval  [0,T].  Since  for  all 
t  €  [0,T],  the  solution  (X0(t),Yo(t),  Z0(t))  remains  in  Og,  the  matrices 

DzMXo(t),Yo(t),Z0(i),0)  and  B(t) 

are  nonsingular  for  all  <  €  [0,T].  | 

Our  aim  will  be  to  determine  when  there  are  solutions  of  (1.1)  that  converge  for  e  — ►  0  to 
a  solution  of  the  reduced  system  (3.1).  For  this  we  introduce  the  following  concept: 

Definition  3.1:  The  system  (1.1a,, b)  is  said  to  degenerate  regularly  on  the  solution  (Xo(t), 
Yo(0)  Z0(t))  of  (3.1),  0  <  t  <  T,  if  a  solution  (x(t,e),y(t,e),z(t,e))  of  (1.1)  exists  on  the  same 
interval  0  <  t  <  T,  which  converges  to  (Xo(<),Vb(f),2’o(0)  as  e  — »  0,  uniformly  on  compact 
subsets  of  0  <  t  <  T. 

The  structure  of  regularly  degenerating  solution  of  (l.la,b)  is  determined  by  replacing 
problem  (l.la,b)  by  two  auxiliary  problems;  the  first  of  these  is  called  the  outer  problem,  and 
the  second  one  the  inner  problem. 
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The  critical  idea  is  here  to  consider  (1.1a)  with  only  an  intial  condition  for  x  but  with  the 
explicit  assumption  that  only  solutions  are  admitted  which  for  e  =  0  reduce  to  a  solution  of 
(3.1).  In  other  words,  we  consider  the  problem: 


X’ =  MX, Y,Z,t), 

*Y'  =  MX, Y,z,e), 

(3.5  a) 

0  =  MX, Y,z,e), 

with  some  initial  condition 

* 

O 

II 

(3.56) 

and  the  limiting  assumption 

*(M>)  = 

x0 «),  Y(t,o)  =  r0(f),  z(t,  o)  =  z0(t) 

(3.5c) 

where  (Xo(t),  Vo(/),  Zo(0)  *s  a  solution  of  (3.1). 

Any  solution  of  (3.5a,b,c)  will  be  called  an  outer  solution.  With  any  such  outer  solution 
we  introduce  in  (1.1a)  the  scaled  variable 


T~tft  (3.6a) 

and  new  dependent  functions 

o(r,  e)  =  t((t,  e)  -  _Y(er,  e) 

/?(T,  =  y(^T,  e)  —  K(ct,  e)  (3.06) 

7( r,  e)  =  z(er,  e)  -  Z( er,  e) 

Let  (i(t,€),y(<,e),z(*,e))  be  a  solution  of  (l.la,h),  then  we  find  that  (a,0, 7)  satisfies  the 
following  DAE,  which  is  called  the  boundary  layer  problem  or  inner  problem: 


^  =  /2(fT,a,/?,7,e), 

0  =  /jfer, «,  (},  7,  e), 

Q(0>  e)  =  t(e)~  C(e),  0(0,  e)  =  77(e)  -  Y( 0,  e),  7(0,  e)  =  C(e)  -  Z( 0,  e), 

where 

/<(*,«,/?,  7,  <0  =  /.(o  +  X(t,e),0  +  y(<,e),7  +  Z(f,«),e) 
-f,(X(t,e),Y(t,e),Z(t,e),e),  i  =  1,2,3. 


(3.7a) 


(3.76) 


To  study  the  asymptotic  behavior  of  solutions  of  (1.1),  we  need  following  assumption 
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Assumption  (II):  The  components  of  the  initial  point  (£(e),t](e),  ((e))  possess  asymptotic 
expansions  : 

oo  oo  oo 

£(«)  ~  Y  »?(e)  ~  Y  W*'  t(e)  ~  Y  c  "♦  °-  (3-8) 

i=0  i=0  »=0 


4.  Asymptotic  expansions  of  outer  solutions 

For  the  analysis  of  the  solutions  of  (3.5a,b)  and  (3.7a,b)  and  their  interrelationship,  asymp¬ 
totic  considerations  are  to  be  used.  We  motivate  here  briefly  the  approach  and  defer  proofs  to 
the  next  subsection.  Suppose  that  the  initial  function  £*(e)  of  (3.5b)  satisfies 

oo 

T(e)  ~  6  +  Y  £*' €'’  as  e  -+  0,  (4.1a) 

i=  1 

and  accordingly  that  any  outer  solution  (X(t,  e),  Y^e),  Z(t,  e))  has  a  formal  asymptotic  ex¬ 
pansion  in  terms  of  e 

N 

X(t,e)  =  YMiV  +  0(eN+1)  as  e  — *  0 

t=0 

N 

Y(t ,  e)  =  Y  y>(0e‘  +  0(eN+1 )  as  e  -»  0  (4.16) 

«=o 
N 

Z(t,e)  =  +  0(eN+1)  as  e  — +  0, 

«=o 

which  is  assumed  to  hold  uniformly  for  0  <  t  <  T,  0  <  e  <  ei  (ei  <  e'). 

Inserting  (4.1b)  into  the  equation  (3.5a),  expanding  the  right  side  functions  at  the  point 
(Jfo(f),yo(t),  Zo(i),0)  and  equating  coefficients  of  equal  powers  of  e,  we  obtain  that  the  first 
term  (Xo(t),  Y0(t),  Zo(0)  must  be  a  solution  of  the  reduced  problem  (3.1a.b)  while  the  higher 
terms  (A"r(t),  Yr(f),  Zr(t )),  r  =  1, ...,  N  in  (4.1b)  must  satisfy  a  linear  DAE  of  the  form 

~  =  fl*(t)Xr  +  fly(t)Yr  +  fu(t)Zr  +  pr(t), 

=  f2*(t)X  r  +  hyWr  +  h,(t)Zr  +  9r(t),  (4.2)r 

0  -  f3x(t)Xr  +  hy(t)Yr  +  f3,(t)Zr  +  rr(i), 

*r(0)  -  C, 


7 


where 


/«(*)  =  ^r/i(^o(O,n(O,Z0(*),0), 

/i»(0  =  £>y/i(A'o(t),y^(0,Zo(0,0),  i  =  1,2,3.  (4.3) 

/»(<)  =  ^*/i(^o(0,n(0,^o(0,o), 

and  the  terms  pr(0»?r(0  and  rr(<)  are  polynomials  in  Ari,Vi,  Zj,  . . A'r_i,  Vr_i,Zr_|  for 
which  the  coefficients  are  higher  derivatives  of  the  functions  at  the  point  (X0(t), 

Vo(t),  0).  The  right  side  £*  of  the  initial  condition  is  the  corresponding  coefficient  in 

the  asymptotic  expansion  of  £*(e)  and  pr,  qr  and  rr,  r  =  0, 1, . . . ,  N,  are  obtained  recursively. 
Therefore,  (pr(t),Qr{t),rr(t))  is  well  defined  on  the  interval  0  <  t  <  T  if  the  previous  terms, 
X] ,  Y\ ,  Z\ XT-\,  Yr-u  Zr-u  are  well  defined  on  [0,T].  A  comparison  of  (4.2)r  shows  that, 
formally,  all  coefficient  functions  satisfy  a  linear  system  of  the  form 

dx 

—  Qn(0x  +  a\i{t)y  +  ai3(0x  + 

0  =  +  a22(t)y  +  a23(t)r  +  62(<),  ^4  4^ 

0  =  fl3i(*)x  4-  o32(t)y  +  033(0*  +  M0> 

x(0)  =  ^, 


for  which 

(a)  a,y,f)i  are  continuous  on  (0,  T]; 

(b)  the  matrix 

Ai(f) 


( °22(0  ^23(0^ 

\  ®32  (0  «33 (0/ 


is  nonsingular  for  t  e  (0,  T). 

For  such  systems  we  obtain  from  Proposition  2.1  the  existence  result: 


Proposition  4.1:  Under  the  assumptions  (a)  and  (b)  the  system  (4.4)  has  exactly  one  solution 
(x(t),y(t),z(t))  defined  on  the  interval  [0,T]. 

Proof:  By  assumption  (b)  we  can  solve  the  second  and  third  equations  in  (4.4)  for  (y,*)  in 
terms  of  x  substitute  into  the  first  equation  in  (4.4).  Thus,  we  obtain  an  initial  value  problem 
for  a  linear  ODE.  Then  applying  the  basic  existence  theory  for  DDEs  (see  Theorem  5.2  in 
[CoLe])  we  know  that  the  system  (4.4)  possesses  a  unique  solution,  g 

Since  the  informal  expansion  procedure  provides  that  pr(f),  gr(f),  rr(f)  are  polynomials 
in  Xi,Yi,Zi,  Ar_i(f),  Fr-i(<),  ZT~\(t)  with  the  coefficients  depending  on  (Ao(t)>  Y0(t), 
Zo(t)),  w<-  cajl  verify  easily  that  the  conditions  (a)  and  (b)  for  the  linear  system  (4.4)  will  hold 
for  all  systems  (4.2)r,  r  =  1 , . . . ,  N.  Then  by  applying  Proposition  2.1  to  these  systems  (4.2)r, 
wc  obtain  the  following  existence  theorem: 
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Proposition  4.2:  If  the  conditions  in  Proposition  3.1  hold  then  each  system  (4.2)r,  r  = 
1, . . . ,  N,  with  any  given  value  of  £*,  has  a  unique  solution  ( Xr(t),Yr(t),  Zr(t ))  defined  on  the 
domain  [0,T]  of  (X0(t),  Y0(t),  Z0(t))- 

Proof:  For  any  r,  1  <  r  <  JV,  the  coefficient  matrix  in  the  system  (4.2)r  , 

(f2y(t)  /a,(t)\ 

is  nonsingular  for  all  t  £  [0,  T]  due  to  the  nonsingularity  of  B(t).  Hence  by  Proposition  4.1,  we 
find  that  the  system  (4.2)r  has  a  unique  solution  defined  on  the  interval  [0,T],  | 

Note  that,  although  the  sequence  of  systems  (4.2)r  is  derived  formally  under  the  hypothesis 
that  their  solutions  (JYr(t),yr(t),  Zr(t)),  r  -  0, 1, . . .  ,N,  are  the  coefficients  of  the  expansion 
series  of  an  outer  solution  (X(t,e),Y(t,e),Z(t,f))  of  (1.1a, b),  these  systems  (4.2)r  themselves 
are  independent  of  the  concept  of  an  outer  solution.  So  far,  we  only  proved  the  existence  and 
uniqueness  of  the  systems  (4.2)r.  Obviously,  this  does  not  mean  the  existence  of  an  outer 
solution.  But  motivated  by  the  procedure  used  in  the  derivation  of  the  systems  (4.2)r,  we  may 
exploit  the  solutions  (Xr(t),  Vr(<), Zr(i)),  r  =  0, 1, ...  ,N,  to  construct  an  outer  solution.  This 
will  be  discussed  in  the  next  chapter. 


5.  Existence  of  outer  solutions 

This  chapter  concerns  the  existence  of  an  outer  solution  of  the  outer  problem  (3.5a, b,c). 
We  cite  the  following  Lemma  which  plays  an  important  role  in  the  study  of  singularly  perturbed 
ODE’s. 

Lemma  5.1:  Let  A(<)  be  an  n  by  n  continuous  matrix  for  t0  <  t  <  t\  and  let  the  real  parts 
of  all  its  eigenvalues  be  less  than  —p  on  to  <  t  <  t\  for  some  p  >  0.  Let  <j>(t ,  s,  e)  be  the 
fundamental  solution  of 

^=A(t)X 
X(a)  = 

on  to  ^  ^  some  s  onto  <  s  <  tj .  Then  there  exists  a  constant  K,  which  is  independent 

of  e,  such  that 

||^(t,s,e)|i  KKe-^-’V2* 

for  to  <s<t<h. 

For  the  proof  see,e.g.  [Le]. 

For  the  theory,  we  require  further  assumptions  about  the  solution  (.Y0(t),  Y0(t),  Z0(t))  of 
the  reduced  system  (3.1a,b). 
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Assumption  (III):  Proposition  3.1  holds  and  for  any  t  G  [0,T]  all  eigenvalues  of  the  matrix 
B(t),  defined  in  (3.3b),  remain  strictly  in  the  left  half  plane. 

In  the  previous  chapter  we  formally  derived  the  systems  (4.2)r  by  assuming  that  an  outer 
solution  (X(t,  e),Y(t,(),  Z(t,e))  has  an  expansion  of  the  form  (4.1b),  and  we  obtained  exis¬ 
tence  and  uniqueness  results  for  solutions  (JVr(t),Kr(t), ZT(t))  of  (2.3)r,  r  =  1,...,AT.  The 
following  theorem  shows  that  these  solutions  (Xr(<),  Yr(t),  Zr(t))  can  be  used  to  construct  an 
approximation  of  an  outer  solution  (X(t, «),  K(t,  e),  Z(t,  ()). 

Theorem  5.1:  Under  Assumptions  (I)  -  (III)  and  for  any  given  £*,  r  —  1, . . . ,  A7,  there  is  a 
constant  f()  >  0  (Vo  G  Jt )  such  that  the  outer  problem  (3.5a,b,c)  has  for  any  f,  0  <  e  <  f0 
a  solution  X  ~  X(t,e),  Y  =  Vr( i ,  e) ,  Z  =  Z(t,t )  G  'D,  defined  on  the  same  interval  [0,T]  as 
(ATo(t),Fo(Oi  ^o(0)>  whirh  satisfies 

N 

X(t,e)-  ^Xi(i)ci  =  0(eN+,)t 

i=0 

N 

Y^(ty  =  0{eN+1),  as  f  — ♦  0  (5.1) 

1=0 

JV 

^(<Ic)-^2ri(0e,'  =  O(eW+,)l 

«=o 

uniformly  for  0  <  t  <  T. 

[Proof]  To  begin  the  proof  we  simplify  the  outer  problem  (3.5)  by  introducing  a  change  of 
variables  defined  by  the  affine  mapping 

T,  :  Rm  x  Rn  x  Rk  -»  Rm  x  Rn  x  Rk  ;  Tt(u,v,w)  =  (X,Y,  Z),  (5.2a) 

N 

x  =  u  +  Yxr(ty, 

r=0 

N 

Y  =  v  +  £rr(f)er  +  Ax(t)u,  (5.26) 

r=0 

N 

Z  —  w  ^  Zr(t)fr  +  Ai(t)ti  +  Bi(t)u. 

r=0 

Here  A\,  A2  and  B\  are  chosen  as 

AAt)  =  u(0"I(/2.(0(/3*(0)"73*(0  -  /2*(0), 

Mt)  =  -(/3,(0)"73,(*Ml(0  -  (/3.(0)"73*(0.  (5-3) 
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where  B(t )  is  defined  in  (3.3b).  Since  the  domain  of  the  functions  /i(:r,  y,  z,  e),  /2(:r,y,z,e) 
and  /3(a:,  y,  z,  e)  is  V  x  J(,  we  require  that  the  new  variables  (u,  u,  w)  remain  in  some  suitable 
domain  such  that  the  range  of  the  mapping  Tt  belongs  to  V  for  all  0  <  t  <  T  and  0  <  e  <  €j, 
where  ei  €  Jt  is  sufficiently  small.  In  order  to  find  such  a  domain  for  Tt  it  is  important  to 
note  that  (Xo(t),Yo(t),  Zo(t))  is  an  interior  point  in  T>  for  any  t  €  [0,T].  In  fact,  since  (5.2b) 
can  be  written  as 


*  =  *„(*)+«+  (e-M'K*1)  «, 

Y  =  Y„(i)  +  v  +  A,(t)u  +  ^  rr(t)er-^  e, 

Z  -  Z0{t)  +  A2(*)u  +  4-  w  +  fy^Zr(t)er~1  J  e. 


(5.4) 


and  (Xr(t),Yr(t),  Zr(t)),  r  =  1, . . . ,  N,  Aj,(<),  A2(t)  and  5i(f )  are  uniformly  bounded  in  t,  0  < 
t  <  T,  it  follows  from  (Jf0(<),lo(0»  Zo(t))  6  int  ( V )  V<  €  [0,T]  that  there  exist  positive 
numbers  €j  (<  e0  G  Jt)  and  S  such  that  Tt(u,v,w)  €  V  for  any  given  t,  t  €  [0,T]  and  for  all 
(u,v,w)  €  5(0, 6)  x  5(0, 6)  x  5(0,  <5).  Accordingly,  for 

V  =  Vu  x  Vv  x  T>w, 

Du  =  5(0,i)cF,  P„  =  5(0,6)cF,  Vw  =  5(0,  tf)  C  Rk, 


we  have  Tt(T>)  C  V  for  all  0  <  t  <  T. 

Substituting  (5.2b)  into  the  first  equation  of  (3.5a)  we  obtain  that 


du  dXAt)  { 

1=0 

N  N 

=  fi(Xo(t)  +  u  4-  Vo(0  +  v  +  Ai(t)u  4  yi(0e*'  (5-5) 

i  i'=i  »=i 

N 

%o(t)  +  w  4  A2(t)u  4  B\{t)v  +  ,  ()■ 

1  i- 1 

I 

I 

; 

!  We  introduce  the  Taylor  expansion  of  f\  at  (X0(<),  Yo(i),  %o (<),  0)  .  For  this  the  abbreviations 
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will  be  used 


fi,(0  = 


=  fi0  +  + 


/x,(t)\ 

2$  •  '■* . * 

V  / 


(5.6) 


where  S^i  is  the  Kronecker  delta.  Then  the  expansion  of  /i(fl)  at  Qo(0  has  the  form 

/.(«) 

N  N 

=  /.(fto(O)  +  E  -0(r,/,(fto(OXE  ft<(0<T  +  Wi(flo(t)M(W  f 

r=l  r'  1  =  1  (57) 

JV  JV  N 

E  H  D<r)/l (fto(0)M(i)W  +  E  n*(0'T  -  (E  +  G(*’ “> f) 

r=2  r'  «=l  i=l 

where  G(t,  u,  u,  w,  t)  is  the  remainder  term.  Obviously,  G(t,  u,  v,  u>,  e)  satisfies  the  following 

Condition  (N) 


Definition  5.1 :  AC 1  -function  F(t,u,v,  re,  f )  :  R  x  Rm  x  Rn  x  Rk  x  R  —*  R*  is  said  to  satisfy 
the  Condition  (N)  if  it  satisfies  the  following  asymptotic  relations: 


0) 

F(<,0,0,0,e)  =  0(fN),  asf-+0; 

uniformly  for  0  <  t  <  T,  where  N  is  a  positive  integer; 

(H) 

DuF{t,u,v,w,e)  =  0(e  +  |u|  +  |u|  -f  |u>|)  \ 

DvF(ty  u,  t>,  w, «)  =  0(e  +  |ti|  +  |t’|  +  |tt'|)  >  as  e,  |u|,  jv|,  Ju>|  — >  0 

DwF(t,u,v,u>,()  =  0(e  +  |u|  +  |v|  +  |u>|)  J 
uniformly  for  0  <  t  <  T. 

It  is  noticed  that  the  Condition  (N)  is  the  same  conditions  as  (1.2a)  and  (1.2b)  in  Hy¬ 
pothesis  (H)  in  [Ya2].  Observe  that 
D/i(fto(OyV 

=  £>,/.(ft„(0)u  +  DyMQMMv  +  yi,(0<i)  +  D*/i(fto(0)(w  +  Mt)u  +  B}(t)  v) 

( D.o  J 

=  (Dt/,(ftu(t))  +  DMSlo{i))Ai(t)  +  DJl(n0{t))A2(t))u+ 

(Dvfi  (ft0(<))  +  ^z/i(fto(f  )Bl  (t  ))v  +  Dx  fx  (fi0(<))u> 
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52  3c<r>/i(S!o(<))(^  «i(')<‘)r 

r=l  r'  i=l 

E  fn  «*,(*))  *‘,+  +‘- 

r=1  '  *1 . fcr-1  \i=!  / 

iV 

=  E*.(0«* +  **(*,«) 


».(*) = e  4j5(r)/i(fi0(<))  e  n  *m*> 


and  R.N(t,  e)  =  0(eAr+1)  is  independent  of  (u,  u,  w).  Since  (X,(t),  Y,(t),  Za(t)),  s  =  0, 1, . . . ,  N 
satisfy 


where  b0(t)  =  /i(fi0(/)),  it  follows  from  (5.5, 6, 7)  that 

^  =  (DMQoit))  +  Dyfi(&o(t))Ai(t)  +  DJ1(Sl0(t))A2(t))u+ 

(Dyfi(Qo(t))  +  +  Z?,/i(£2o(<))»  + 

+ E  +f:my)r  -  (^(tyn  (5-9) 

r=2  r'  t=l  .=  1 

+  <?(<,  U,  V,  W,  e) 

==  Ci(<)u  +  Id(t)i>  +  £?i(<)«»  +  Fi(i,  U)  v,  w,e) 


C,(t)  =  Dx/i(«o(0)  +  D,/i(no(i)Mi(i)  +  ^«/i(no(0M2(<), 
il(0  =  DyMSlAt))  +  D,/,(tto(<)-Bi(0. 

Exit)  =  2?,/,(fto(0), 

/V  N  N 

=  £  ^D<r>/J(no(*))[M(<)W  +  Efi«(OOr  -  (E^KD 


Obviously 


+  Rw(t,  e)  4-  G(f,  w,  v,  w,  e). 


RN(t,  e) 


and 


N  .  N  N 

E  ^(r)/i(n0(<))[(„4(*)^  +  E  n.(OOr  -  (EnW)r] 

r=2  r'  i=l  «=1 

satisfy  Condition  (N).  Since  also  G  satisfies  Condition  (N),  this  implies  that  Fy  satisfies 
Condition  (N)  as  well. 

Now  substitute  (5.2b)  into  the  second  and  third  equations  of  (3.5a)  and  expand  /2,  f3 
at  (X0(<),  Vo(t),  Zo(t),  0).  Then,  because  (A\(i)>  Vj(<),  Zi(t)),  i  =  0,1,...,  Ar,  are  solutions 
of  (4.2 )i,  we  obtain  two  further  equations  which,  with  (5.9)  form  the  the  following  system 
equivalent  to  (3.5)  under  the  mapping  %\ 

du 

—  =  C,(t)u  +  £i(<)t>  +  &i(0u’  +  Fi(t,u,  v,te,e) 

e—  =  B(t)v  +  E2(t)w  +  u,  t>,  7 r,  e) 

0=  E3(i)iv  +  F3(t,u,t;,te,e)  ^  ^ 

N 

«( 0,0  =  r(«)  -  E^r  =  MO  =  0(eN+1) 

r=0 


Here 


E2(t)  =  Dtf2(Q0(t)),  E3(t)  =  Df/,(Q0(<), 

and  F2  and  ,F3  satisfy  Condition  (N).  Note  that  in  the  derivation  of  the  last  two  equations 
of  (5.11),  we  used  the  fact  that 


DzMQ0(i))  +  D,/2(flo(OM,(0  +  DJ2((l0(t))A2(t)  =  0, 

Dzf3(nQ(t))  +  D,Mn0(t))A,(t)  +  £>,/3(fio(OM-2(0  =  o,  (5.12) 

DMn0(t))  +  DtM^{t))Bi(t)  =  o, 

and 

Dyf2(n0(t))  +  Dtf2(Cl0(t))Bi(t)  =  i?(<) 

when?  B(t.)  is  defined  in  (3. 3b). With  the  notations  (5.6)  this  is  equivalent  with 

Df2(n0{t))A(t)  =  (o,B{t),Dj2(n0mDj2(n0(t))), 

^/3(fio(OM(0  =  (o,o,pz/3(n0(#)),A/3(fio(0))- 

In  other  words,  from  the  defining  relation  (5.3)  we  find  that  Ai(t),  A2(t)  and  Bi{t)  solve  (5.12). 

If  we  can  show  that  there  exists  a  solution  (u(f,  e),  v(t,  e),  tv(t,e))  6  P  of  the  system  (5.11) 
which  satisfies 


u(t,e)  =  0(eN+')y  v(t,e)  =  0(eN+i ),  «.(f,f)  =  0(eN+1), 


(5.13) 
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uniformly  for  0  <  t  <  T,  then  Theorem  5.1  is  proved.  But  note  that  we  need  only  the  existence 
of  a  solution  of  (3.5)  without  requiring  its  uniqueness.  Hence  instead  of  (5.11)  we  consider  the 
constrained  system  of  integral  equations 

u(t,  e)  =  $(t)(Me)  +  J'  *-Hs)(Li(s)v(s,  e)  +  F,(s)u>(s,  e) 
v(t,e)  = 

w(t,e)  =  -(F3(i))_1F3(i,u(t,e),v(i,e),tu(t,e),  e) 
where  $(t)  satisfies 

2M  =  C,  «*(«).  0<t<T 

$(0)  =  I 

while  ^(tjS,  e)  is  the  solution  of  the  following  system: 

1 

-T  =  ±B(t)V,  0  <  a  <  t  <  T 

at  e 

*(t,s,e)\t=a  =  I. 

Note  that  the  system  (5.14)  is  not  equivalent  with  the  system  (5.11).  But,  obviously,  if  (u(t,  e), 
v(t,  e),  w(t,  e))  G  V  solves  (5.14),  then  it  solves  (5.11)  as  well.  Thus  if  we  can  prove  that  (5.14) 
admits  a  solution  which  satisfies  the  asymptotic  relation  (5.13),  then  we  are  done.  We  will  use 
the  theorems  in  [Ya2]  to  prove  this  existence. 

By  Assumption  (III)  there  exists  a  positive  number  \i  such  that 

«(*(*(*)))  < 

for  all  eigenvalues  of  B(t),  0  <  t  <  T.  Then  Lemma  5.1  ensures  the  existence  of  a  constant  K, 
which  is  independent  of  e,  and  such  that 

WM,e)|  0<s<t<T. 

With  the  notations 

0i(*,e)  =  *(<)M«).  02(t,e)  =  0,  (s)Li(s), 

H^t,s,€)  =  ^(t)^~1(s)Ei(3),  K(t,s,€)  =  V(t,s,e)E2(s),  (5.15) 

F]  =  $(*)$“x(3)jFi,  F2  =  F2,  F3  -  -E~\t)F3, 

(5.14)  can  be  written  in  the  form  of  constrained  systems  of  integral  equations  as  the  one 
in  [Ya2],  and  evidently  the  Hypothesis  (H)  in  [Ya2]  are  satisfied.  Therefore  by  the  results  of 


+  Fi(s,u(s,  e),  u(s,c),u;(s,e),  e))ds) 
f  e\E2(s)w(s,  e)  +  F2(s,  u(s,  e),  v(s,  e),  u>(.s,  e),  e))ds 
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Theorem  1,  2  in  [Ya2],  (4.13)  has  a  unique  solution  (u(t,  e),  v(t,  e),  u?(t,  <?))  e  T>,  for  all  t  €  [0,X] 
and  0  <  f  <  (2  (  <  fi)  that  satisfies  (4.12).  Therefore,  under  the  transformation  (4.1),  (A'(f,  e), 
Y(t,e),  Z(t,e))  remains  in  T>  for  all  t  €  (0,T]  and  0  <  e  <  e2  and  is  an  outer  solution  which 
satisfies  the  asymptotic  relations  (4.1b).  This  completes  the  proof  of  Theorem  5.1.  | 


6.  Asymptotic  expansion  of  inner  solutions 

6.1.  Formal  derivation  of  asymptotic  expansion  of  inner  solutions.  Now  to  proceed 
to  the  inner  system,  by  using  Assumption  II,  for  the  initial  conditions  o(0,  e),  0(0,  e),  7(0,  c) 
in  (3.7a)  we  must  have 

N 

a(o,  C)  =  f(f)  -  no  =  £(6  -  av  + 

«=i 

N 

0(0, 0  =  vie)  -  r(0, 0  =  J>.  -  Y.( 0)K  +  0(eA,,+1 ),  (6.1) 

1=0 

7(0. «)  =  C(<)  -  Z( 0,0  =  £«.  -  Zi(0))e'  +  0(e"+1). 

»=0 

Then,  formally,  we  expect  the  solutions  of  (3.7a)  to  satisfy  asymptotic  relations  of  the  form 

N 

«(r,0  =  £oi(r)f''  +  0(fN+1) 

i-0 

N 

0(r,e)  =  ^0,(rK  +  O(£N+1)  (C.2) 

t=0 

N 

7(r,0  =  £7.(r)e,  +  O(fN+1) 

«=o 


which  should  hold  uniformly  for  0  <  r  <  T/e. 

We  will  substitute  (6.2)  into  (3.7a),  expand  right  sides  at  (Xo(D)  +  oro('r).  l’o(0)  +  0o(r)> 
Zo(0)  +  7o(f))  and  collect  terms  with  equal  powers  of  c  to  obtain  the  equations  which  (o,(r), 
0i(r) ,  7i(r))  must  satisfy.  Since  the  expansions  of  the  right  sides  in  (3.7a)  involve  very  tedious 
derivations,  we  set  for  abbreviation 


(o(r,f)  +  X(er,  e)  \ 
flir^e) +  Y(er,e) 
7(r,e)  +  Z(eT,e) 


,r(r,C)  = 


fi(r,e)  = 
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A'(er,e) 

Yicr.e) 

Z(eT,e) 


Then 


/,(er,  a(r,  e), 0(t,  e),  j(t,  e),  e)  =  fi{Cl(T,  e))  -  /«(r(r,  e)),  i  =  1, 2, 3. 

For  simplicity,  in  the  following,  the  subscript  i  will  be  omitted. 

For  later  use  we  state  the  following  two  identities: 


X>;e‘+’’  =  £<!>*-<)«*  +°('N+") 


M*)»  E***1  -  E 

\fc=i  /  Jblt...,k„=i»=i 

From  (6.4)  it  follows  that  in  the  polynomial  d„h„(e )  in  e  the  coefficient  of  er  is 


Y,dn  Y  IK.  ifO  <r<N 


while  for  r  exceeding  N  it  is 


Ea»  E  IK 


Since  (X(t,  e),  T(t,e),  Z(t,e))  and  (a(r,  e),  /3(r,  e),  7(r,  e))  satisfy  the  asymptotic  relations 
(4.1a, b)  and  (6.2),  respectively,  we  obtain,  with  (6.3),  that 
X(er,  e)  +  a(r,e) 

=  ^o(^)  +  a0(r)  +  J(Ti(tr)  +  ^(r))^  +  0(eN+!) 

tl  (6-7) 

N  k  y(*-0/nVr*-i 

_  v  /n\  i  i  \  i  Ai  vuJr  i 


:  ^o(O)  +  “o(t)  +  ^(ajfc(r)  +  Y  '  (l  _  t\i - )€*  +  ^(eW+1)> 


Jfc=l  i=0 


and  similarly  that 


X(er,e)  =  X„(0)  +  £(£  +  0(eN+1).  (6.8) 

With  (6.7),  (6.8)  and  the  analogous  expansions  for  Y(er,  e)  +  a(r,  e),  F(er,e),  Z(er,  e)  +  a(r,  e) 
and  Z(er,e),  we  find  that 

N 

«(t,  e)  =  SI  o(t)  +  5]  nt(r)<*  +  (eN+l), 


r(r, «)  =  r0  +  r,(j-)e*  +  (eA 
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where 


/Xo(0)  +  o0(t)\ 
on/Tx_  W)  +  /?0(r) 
n°(r)"  ^Zo(0)  +  7o(r)  y  nk{T)~ 


AW  +  Ef.. 


7*(t)  +  E,=o  ‘  • 


0  (t-i)! 

(*:->)! 


r»=  ^(O)j'  r‘(r)  = 


<11 

1  L,=0  (i-.)! -  ' 

yL(*“°(0)r*-’ 


(6.10) 


V 


=0  (*-«>• 

Zj*“°(0)r* 

2^.=0  (i  — i)! 

6*1 


(6.11) 


(6.12) 


and  6^1  is  the  Ivronecker  delta.  It  follows  from  (6.10)  and  (6.11)  that 

;(er,a(r,€),0(T,f),7(r,e),e)  =  /(fi(r,e))  -  /(r(r,e)) 

N  N 

=  /(flo(r)  +  Y,  n^r)ek  +  0(tN+'))  -  f(T0  +  £  r*(r)efc  +  <V'+1 )) 

*=]  *=i 

i=l  \Jt=l  / 

i=l  -7'  \*=1  / 

By  collecting  equal  powers  of  e  in  (6.12),  with  (6.6)  we  find  that  the  coefficient  cr(r)  of  er, 
1  <  r  <  TV,  is: 

«,<r)  =  £^f^  £  ri«,M-£^>  £  nr,.,,) 

>=  1  *,+  +*>->■  »=1  J  =  1  ‘l+  +*,='J=1 


=  D/(fio(r))nr(r)  +  ^£V<^!l»  £ 

j=2  3'  ‘1+  '+*>=’■  a  =  l 

1<‘>  <>• 

-£^  £  flr,». 

j~  1  *1  +  +*>-■■  S=1 


(6.13) 
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Since,  by  (6.10), 

^/(fto(r))ftr(r)  =  nxf(Q0(r))ar(r)  +  D,f  (Sio(r))fir(r)  +  DJ(fl0(r))7r(r) 

+  D,mo(T))±  +  Wr))t  (,14) 

+  D,m{T))  £  +  D.m „M)<V, , 

it  follows  from  (6.13)  that 

cr(r)  =  JD3./(J20(T))or(r)  +  £>j,/(ft0(r))/?r(r)  +  DJ(SlQ(T))jr(T)  +  Pr(r)  (6.15) 


=  D„/(ft,(r))£ 


+  fl,/(n„(r))  £  +  D,/(fio(r))«rl 

i=0  ' 

E  E  ^rt( 
j—2  *1+  +*;=>•  s=l  i- 1  *!+■  +*j=r  *=1 

x^r-,Voirfc”' 

=  (Dtf({lo(r ))  -  D,/(r„))tf,,  +  (D,/((i„(r))  -  D,/(r„))  £  - 

+  (D,f(a0(r))  -  D,m))  £  ’ 

+  (D./(fio(r))  -  D,/(r0))  £ 


+f£m mi  £  £  nr*.M 

j  =  2  »1+  +kj=r  3=1  i  =  2  *,+  +kj=rS=l 

l<*j<r  »<*><•■ 

which  is  a  polynomial  in  A*-i,7r- 1  with  the  coefficients  depending  on 

£>7(fto(r)),  j  =  l,...,r  and  £P/(r0),  (xf(0)t  V/°(0),  zJO(0)),  i,j  =  0,1, The 
constant  term  in  the  polynomial  Pr(r)  is 

E(oi/(!joW)-£’,«r»))/j!  e  nr*.M 


f'  C'/(r0) 


=  Y,  (®WW)  -  DV(r.))  PrjW 
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where 


Mr)  =  7i  E  (6.16c) 

*l4-+»j.r  *  =  J 
l<ki<T 

Obviously,  by  (6.11),  the  degree  of  the  polynomial  prj(r)  in  r  is  <  r.  In  other  words,  the 
constant  term  in  the  polynomial  PT  is  the  sum  of  certain  terms,  each  of  which  is  a  product  of 
the  factor  of  the  difference  between  the  derivatives  of  /  at  ft0(r)  and  T0  and  the  polynomial 
p,j(r)  in  r  with  degree  <  r. 

Now  substituting  (6.2)  into  (3.7a),  with  (6.15)  wc  obtain  that 

+<*»+■> 

r=:  0 

=  f(/.(fi(v))-/>(r(r,e))) 

=  <{/i(fto(r))  -  /i(r0(r))  +  fc'r»(r)0  +  0(eN+1), 


N 


Y 


d0r{T) 

dr 


(r  +  0(eN+1 


) 


N 


=  /3(0o(r))  -  /2(r„)  +  Y  4a)(0er  +  0(eN+1 ), 

r=l 

N 

o  =  /3(fio (r))  -  /3(r„)  +  £  43>(r)er  +  0(e*+I ), 

r=1 

where  4'*(t)  is  defined  in  (6.15)  with  /  replaced  by  /;,  i  =  1,2,3.  Since 


/a(r0)  =  /2(xo(0),ro(0),z0(0),0)  =  o, 

/3(r0)  =  /3(Ar0(o).r0(o),z0(o),o)  =  o, 

we  find,  by  collecting  terms  with  equal  powers  of  e,  that  a0(r)  =  0,  because  =  0,  ft0(0)  =  0, 
and  (  since  Ao(0)  =  £0,  ) 


=  /a(6,  Y0(0)  +  0o,  Zo(0)  +  7o,  0), 

0  =  /3(6.>o(0)  +  0o,Zo(O)  +  7o,0),  (617)o 

A»(0)  =  7o-  V'o(0),7o(0)  =  Co  -  -Zo(O), 
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and,  generally, 


=  Dxf2(Sl0(T))ar  +  D,/a(no(r))j9r  +  P*/2(fto((r))7r  +  ^,r(r),  (6.17)r 

0  =  D,/s(no(r))ar  +  Dyf3(n0(r))Pr  +  DJ3(n0(r))7r  +  P7,r(r), 

ar(0)  =  Cr-C^r(O)  =  T}r  ~  Vr(0),7r(0)  =  Cr  ~  Zr(0) 

for  r  =  1, . . . ,  N,  where  Pa,r(T )  =  cr-i(r)»  anc*  Pp,r(T)i  ^7 ,r(T)  are  defined  by  (6.16a,b,c)  with 
/  replaced  by  /2  and  f3 ,  respectively. 

The  systems  (6.17)r,  r  =  1, . . . ,  N,  are  all  linear  systems  in  cvr,  /3r,  7r.  In  particularly,  the 
right  side  of  the  first  equation  of  (6.17)r  is  independent  of  ar,  /3r,  7r,  which  means  that  or(r ) 
can  be  obtained  easily  from  the  first  equation  of  (6.17)r.  The  terms  Pa,r(r),  Pp,r(T)’  PfAT)  *n 
(6.17)r  are  polynomials  in  ai,/?i,7i,...,ar-i,/3r- i,7r- 1*  Moreover,  the  coefficients  of  PQr(T) 
depend  on  r,  A)(r),7o(T),  and  (-Xo(0),yb(0),Zo(0))  and  the  derivatives  of  (-X'j(f),  Yi(t),  Zi(t)) 
at  t  —  Ofor  t  =  0, 1, . . .  ,r  —  1,  while  the  coefficients  of  PptT(r),  P7)r(r)  depend  on  t,  /?o(t),7o(t), 
and  (-Xo(O), Fo(0),  Zo{0))  and  the  derivatives  of  ( Xtft ),  Yj(<),  Z,(t))  at  t  =  0  for  i  =  0, 1, . . .  ,r. 
Hence  for  r  =  1, . . . ,  N,  P0,r(r ),  P^r),  P7lr(r),  are  known  recursively,  if  (or_i,  /3r- j,  7r_i) 
can  be  obtained  recursively  for  known  (JVj(f),  ^(f),  Zi(t)),  i  =  0, 1, . . . ,  N  .  The  constant  terms 
of  these  polynomials  are  the  sum  of  certain  terms,  each  of  which  is  a  product  of  the  factor  of 
the  difference  between  the  derivatives  of  /,•  at  (£o>  To(0)  +  0o (r),  Zo(0)  +  7o(r),0)  and  (£0, 
y0(0),  Zo(0),O)  and  some  polynomial  q(r)  in  r  with  degree  <  r. 

6.2.  Properties  of  the  expansion  of  inner  solutions.  With  the  outer  solution 
(X(t,e),y(f,€),  Z(t,e))  of  Theorem  5.1  we  can  apply  the  change  of  variables  (3.6a, b)  and 
construct  the  boundary  layer  problem  (3.7a,b).  For  the  study  of  the  inner  problem  (3.7a,b)  we 
impose  the  following  assumption  on  the  solutions  of  (6.17)o: 

Assumption  (IV): 

(i)  The  initial  value  problem 

=  m 0.  W)  +  /?,  Z0(0)  +  7, 0), 

0  =  /»(£ 0,  Yo(0)  +  p,  Zo(0)  +  7, 0),  ((US) 

P(0)  =  bt  7(0)  =  c, 

corresponding  to  (6.17)o  has  for 


b  —  T)o  —  Fo(O)  and  c  =  £0  -  Z0(0) 
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a  solution  (/9o(i")>  7o(t"))  defined  on  [0,  oo)  such  that 

fio(r)  =  (£o,*o(0)  +  Po(t),Z0(0)  +  7(l(r))  €  V  for  all  r  >  0,  and 


lim  (fio(r)iTo(’'))  =  (0,0). 

r— »oo 

(ii)  The  matrices 
DJ3(Q  0(r), 

C(r)  =  Dyh(^T))  -  Dr/2(n„(r))(Dl/3(fio(T))r1  D,M<1 0(r)) 
are  nonsingular  for  r  >  0. 

Under  the  assumptions  (I)  -  (IV)  and  with  properly  chosen  £*,  r  =  l,...,Ar,  we  shall 
prove  in  Proposition  6.1  that  the  system  (6. 1 7)r  has  a  unique  solution  (ar(r), /?,.(t),  7r(r)) 
defined  on  [0,  oo),  r  =  1, . . . ,  N,  and  that  this  solution  decays  to  zero  exponentially  as  r  — *  oo. 

The  following  result  provides  sufficient  conditions  under  which  Assumption  (IV)  can  be 
derived  from  Assumptions  (I)  -  (III): 

Lemma  6.1:  If  for  the  initial  conditions  (£(e),  q(e),  C(c))  (II)  the  point  (£o,f?o>Co)  is  suffi¬ 
ciently  close  to  the  point  (£o,Vo(0),Zo(0)),  then  Assumption  (IV)  is  a  consequence  of  Assump¬ 
tions  (I)  -  (III). 

The  principal  part  of  Lemma  6.1  is  a  direct  consequence  of  Lemma  6.3  below  which  in  turn 
can  be  proved  by  means  of  the  following  result: 

Lemma  6.2:  Let 

rfH'4  *  +  /«,*) 

where  A  E  Rn  x  n  is  a  real  (  constant  )  matrix  for  which  all  eigenvalues  have  negative  real  parts. 
Let  f  be  real  continuous  for  (x,  t)  6  i?(0, 6)  x  R+ ,  where  =  {t  |  t  >  0  },  and  8  >  0  is  a 
small  number,  and  let 

f(t,x)~o( |x|)  as  |x|  — >  0 

uniformly  in  t,  t  >  0.  Then  the  identically  zero  solution  is  asymptotically  stable. 

For  the  proof  see,  e.g.,  [CoLe], 

Lemma  6.3:  Let  Assumptions  (I)  -  (III)  hold  and  consider  the  initial  value  problem  (6.18) 
where  the  initial  point  (b,c)  satisfies  the  constraint  condition  /3(£ o>  V0(0)  +  h,  Zo(0)  +  c, 0)  = 
0,  and  (Xo(t),Y0(t),  Z0(t))  is  the  solution  arising  in  Assumption  (III).  Then  for  sufficiently 
small  |fc|,  |e|,  the  DAE  (6.18)  has  exactly  one  solution  (/?(r),7(r))  defined  on  [0,  oo),  for  which 
lo(0)  +  P(t)  e  Vy,  Z0( 0)  +  7(7)  6  T>2  and 

UmOS(r),7(r))  =  (0,0). 
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Proof:  By  Assumption  (I)  the  matrix  Dzfs(£o,  lo(0),  Zq(0),0)  is  nonsingular,  (Yo(0),  ^o(O)) 
is  a  point  of  the  open  set  Dy  x  Dz,  and  we  have  /3(fo,^o(0),  Zo(0),  0)  =  0.  By  the  implicit 
function  theorem  there  are  two  positive  scalars  8p  and  8y  such  that  Yo(0)  +  B(0, 8$)  C  T>y, 
Zo(0)  +  B(0,87)  C  T>z,  and  that  for  any  €  £(0,  fy)  there  exists  exactly  one  7  =  <f>(0)  G 
£(0,  <57)  for  which  /3(£o,lo(0)  +  0,ZO( 0)  +  7,0)  =  0.  Since  the  functions  /2  and  /3  were 
assumed  to  be  sufficiently  smooth,  we  may  expand  them  at  the  point  (£0,  Po(0),  Z0(0),  0)  and 
write  (6.18)  as 

=  Dyf2(£o,Y0(0),Z0(0),0)f3 


+  Dyh(So,Yo(0),  Zo(0),0)7  +  W,7)  (6.19) 

O  =  Dyh(^,Yo(O),Zo(O),O)0 
+  DMZo,  y0(0),  Zo(0 ),  0)7  +  G(P,  7) 

where  F(/l,7),  G(/J,7)  =  0( |/3|2  +  I7I2)  as  j3  — >  0  and  7  — ►  0.  Prom  the  second  equation  of 
(6.19),  it  follows  that 

-r  =  -(D.A(&,n(0),Zo(0),0))-1fl,/3«o,Ko(0),Z„(0),0)^ 

+  0(L8|2  +  ItI2) 

which,  with  <$(0)  =  0,  implies  that 

7  =  <M?)  =  <m  as  »  0.  (6.21) 

Substituting  (6.20)  into  (6.19)  and  using  the  asymptotic  relation  (6.21)  we  find  that 

^  =  £/?  +  £(£).  (6.22) 


R(0)  =F(^,7)-D9/2(eo,Vo(0),Zo(0),0)(Di/3(^,n(0),Z0(0),0))-IG(/?,7) 

=om2  +  h\2)  =  omn  as/?-*o 

and  B  =  £( 0)  is  the  matrix  defined  in  Proposition  3.1  for  which  by  Assumption  (III)  all 
eigenvalues  have  negative  real  parts.  Thus  by  Lemma  6.2,  0  6  £(0, 8$)  is  an  asymptotically 
stable  solution  of  (6.22),  namely,  for  any  0  <  e  <  8p,  there  exists  8  <  6p  such  that 

10(0)1  <  8  (  or  |6|  <  8) 

implies  that  the  solution  /?(r)  of  (6.22)  satisfies 


!(t)|  <  e  for  all  r  >  0 
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and 

lim  /3(t)  ~  0. 

Because  /?(r)  remains  in  the  hall  8(0,6/})  for  all  r  >  0,  if  follows  from  the  continuity  of  <f>  and 
r^(0)  =  0,  that  7 (r)  :=  <f>(ft(r))  €  8(0,6 for  all  r  >  0  and 

lim  7(7-)  =  0. 

This  completes  the  proof  of  Lemma  6.3.  | 

Proof  of  Lemma  6.1:  We  apply  Lemma  6.3  to  (6.17)o  to  complete  the  proof  of  Lemma 
6.1,  Lemma  6.3  ensures  that  when  (r/0,(o)  is  sufficiently  close  to  the  point  (Yb(0)t  Zo(0))  then 
|Ai(0)|,  |7o(0)|  are  sufficiently  small  and  (6.17)0  has  exactly  one  solution  (/?0(r),  70(7"))  defined 
on  (0,oo)  which  satisfies 

rlim(A,(r),  7o(r))  =  (0,0) 

and 

(V'o(O)  +  0o(t),Zo(O)  +  70(r))  eVyxV„  Vr  6  (0,oo). 

It  only  remains  to  prove  that  the  condition  (ii)  of  Assumption  (IV)  holds.  For  this  we  can 
choose  sufficiently  small  8$  and  6 7  such  that  for  all  (/?, 7)  €  B(O,80)  x  8(0,6-,)  the  matrices 

P*/3«o,ro(0)  +  /?,Zo(0)  +  7,0), 

Dyf7-  DJ2(DJ3r1Dyf3^o,Yo(0)  +  0,Zo(O)  +  lyO), 

arc  nonsingular.  This  is  guaranteed  by  the  continuity  of  these  two  matrices  and  their  nonsin¬ 
gularity  at  (/?,7)  =  (0,  0).  Then,  for  f  =  8p  there  exists  a  8  >  0  such  that 

l/?o(0)|  <  8,  1 70 (0)|  <  6, 

implies  that  (/?o(T)i  7a(r))  €  8(0,6)  x  8(0,6)  for  all  r  >  0.  Thus  on  the  solutions  0q(t),  7o(r) 
the  matrices  (6.23)  are  indeed  nonsingular.  This  completes  the  proof  of  Lemma  6.1.  | 


Lemma  6.4:  Under  Assumptions  (I)  -  (III),  if  the  solution  (^o(r)»7o(r))  °f  (6.1 7)0  satisfies 
the  Assumption  (IV),  then  there  exists  a  positive  o  such  that 

fio(r)  =  0(e"r)>  7o(r)  =  ()(*-”)  (6.24) 


as  t  — ►  00. 

Proof:  Since  /?0(r)  —>  0,  70(r)  — *  0  as  r  — ►  00,  for  any  given  e  >  0,  we  may  choose  a 

sufficiently  large  ru  >  0  such  that  /30(r)  €  .8(0,6/})  and  70(7")  €  8(0, 6-,)  for  all  r  >  r0,  and 

WdM)I  <  * 
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where  the  function  R  is  defined  in  (6.22).  The  equation  (6.22)  can  be  written  as  the  integral 
equation 

A»(r)  =  e"B^o(ro)  +  f  e~B^R(0o(s))ds.  (C.25) 

Jt0 

Since  the  real  parts  of  all  eigenvalues  of  B  are  negative,  it  follows  from  Lemma  5.1  with 
A(<)  =  B  and  e  =  1,  that  there  exist  positive  numbers  I\  and  p  such  that 

j|e_flr|[  <  Ke-*r,  for  all  r  >  0.  (6.26) 

Note  that  here  we  use  the  fact  that  the  fundamental  solution  s)  of  the  equation 

is  (j>(r,s)  =  e-(r-4)fl5  for  0  <  s  <  t  <  oo.  Hence  (6.25)  and  (6.26)  together  show  that 

|A»(r)|  <  c-»r\fofa)\  +  'K  e-«r-‘)Ws)\d3 

Jto 

or  r 

eHA>(r)|  <  l^o('ro)!  +  eK  /  e*‘\fo{s)\d3 
Jto 

By  Gronwall’s  inequality  it  follows  that 

\(3o(t)\  <  \Mro)\e-(Kr<>e-^-tK)T  (6.27) 

and  with  e  =  ^  and  cr  -  p/2  we  obtain  from  (6.27)  that 

/?0(r)  =  0(e~<rr )  as  r  — ♦  oo 

which,  together  with  (6.21),  completes  the  proof  of  Lemma  6.4.  g 


Lemma  6.4  shows  that  the  first  term  of  an  inner  solution  of  (3.7a,b)  is  negligible  out¬ 
side  the  boundary  layer,  which  also  is  a  property  of  any  inner  solution  of  (3.7a, b).  Since 
(ar(r),  ($r ( r ) ,  7r(r)),  r  =  0, 1, . . . ,  JV,  were  derived  as  coefficients  of  an  inner  solution  of  (3.7a, b), 
we  expect  that  all  terms  (o;r(r),J3r(r),7r(r)),  r  =  0,1,..., IV,  possess  this  property,  namely 
that  ,  for  r  =  1, . . . ,  JV, 

JKm^  ar(r),  /?r(r ),  7r(r))  =  (0, 0, 0). 

This  will  be  confirmed  in  the  following  Proposition  6.1  where  also  the  existence  and  uniqueness 
of  solutions  (ar(7')?/?r(T),  7r('r))  of  the  systems  (6.17)r,  r  =  1, . . . ,  JV,  is  discussed. 


By  (6.17)r  all  coefficient  functions  («r,  /?r,  7r)  satisfy  a  linear  system  of  the  form 

da  _  .  . 

*=P,(T)’ 

^  =  -42l(r)a  +  A2i(r)fi  +  /lj3(rh  +  C2(r),  (6.28) 

0  -  431(t )ft  +  432(t)/?  +  A33(t)i  +  P3(t), 
a(0)  =  «*,  /?(0)  =  /?*,  7(°)  =  7*, 


for  which 

(i)  Aij,Pj  arc  of  class  C1  on  (0,oo)  and  -A33(t)  is  nonsingular  for  all  r  >  0; 

(ii)  Aij  are  bounded  uniformly  in  r  €  [0,oo)  and  the  limit 

lim  (A22(r)- ^(r^VMsafr))  =  B0 

T  — »00 

exists  ; 

(iii)  the  limiting  matrix  Bo  has  all  eigenvalues  remaining  strictly  in  the  left  half  plane. 
And  the  initial  point  (ft*,  /?’,  7*)  satisfies  the  compatibility  condition 


•431  (0)o*  +  432(0)/?*  +  433(0)7*  +  P 3(0)  —  0. 


For  the  system  (6.28)  we  obtain  the  following  result. 

Lemma  6.5:  Under  the  conditions  (i),  (ii)  and  (iii)  the  system  (6.28)  has  exactly  one  solution 
(<v(r),  /?(r),  7(r))  defined  on  the  interval  [0,  00).  Moreover,  the  asymptotic  relations 


Pi(r)  =  0(e~<Tr),  as  t  -*  00,  i  =  1,2,3 

hold  with  some  er  >  0.  Then  x 

o*  =  -  /  Pj(s)ds  <  00 

Jo 

and  with  this  ft*  there  exists  a  solution  (o(r),  /?(r),  7(r))  of  (6.28)  that  satisfies 
ft(r)  =  0(e~<7' r),  (3(t)  =  0(e'ffir)i  7(0  =  0(*"',r),  as  r  -  00. 


(6.29) 


(6.30) 


with  some  0  <  oq  <  <7  . 

Proof:  From  the  first  equation  of  (6.28)  wc  obtain 

ft(r)  =  a*  +  /  P,(s)ds.  (6.31) 

Jo 

Hence  since,  by  assumption  (i),  43s(r)  is  invertible,  it  follows  from  the  third  equation  of  (6.28) 
that 

7  =  — 4J‘3,(t)  (A3i(t)o  +  432(r)/?  +  P3(t)).  (6.32) 
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Substituting  (6.31)  and  (6.32)  into  the  second  equation  of  (6.28)  we  obtain 

%  =  M  t)0  +  Kt) 


M?)  -  A22 M  -  A23(t)A331(t)A32(t), 

6(r)  =  P2(t)  -  A23(t)A;3\t)P3(t)  +  (A21(t)  ~  ^(r^V^i^Mr). 

Since  (6.33a)  is  a  linear  ODE  with  the  coefficients  defined  on  [0,  oo)  it  follows  from  the  basic 
existence  theorem  for  initial  value  problems  of  linear  ODE’s  that  with  the  initial  condition 
/?(0)  =  /?*  the  equation  (6.33a)  has  a  unique  solution  fl  =  /3(r)  on  [0,  oo).  Inserting  fi  =  /3(t) 
into  (6.32)  we  obtain  the  component  7  =  7(7-)  of  the  solution  for  (6.28).  This  proves  the  first 
part  of  Lemma  6.5. 

For  the  rest  of  this  lemma,  note  that  because  of 


Pj(t)  =  0(e  trT),  as  r  ->  oo, 


the  integral 


r  Pi(s)ds 
Jo 


exists.  Then  by  choosing  a*  as  in  (6.29)  we  find  that  the  solution  <y(r)  has  the  form 


°(T)  /  Pl^S 


which  implies  that  a(r)  satisfies  the  asymptotic  relation  (6.30).  Since  Aij,  i  =  2, 3,  j  —  1,2,3, 
axe  bounded  uniformly  for  r  >  0,  and  P2(t),  P3(t),  a(r)  satisfy  the  asymptotic  relation 

P2(r)  =  0(e”‘rr),  P3(r)  =  0(6"").  *00  =  0{e^r),  asr^oo 

it  follows  that  b(r)  in  (6.33b)  satisfies  the  same  asymptotic  relation,  namely 

b(r)  =  0(e~ar).  (6.34) 

Since  by  conditions  (ii)  and  (iii),  limr_oo  A(t)  —  Bo  and  all  eigenvalues  of  Bo  remain  strictly 
in  the  left  half  plane,  there  exist  /*'  >  0  and  r0  >  0  such  that  for  any  eigenvalue  A(  A(r))  of 
A{t)  the  inequality 

3?(A(A(r)))  <-/*'<  0,  V  r  >  r0 
holds.  Let  4>(r,  s)  be  a  fundamental  solution  of  the  equation: 


$(t,  s)lr=,  =  I,  To  <  S  <  T. 


27 


From  Lemma  5.1  with  e  =  1  we  obtain  the  estimate 


Then  for  the  solution  /?(r)  of  (6.28),  it  follows  that 

/?(r)  =  $(r, 0)/?*  +  J  4>(r,  s)i>(s)cfs 

=  0(c~'l'r/2)  +  j  0(e-'‘'(r-,)/2)0(e-ffi’)^ 

=  0(c"'/r/2)  +  0(c~aT) 

=  0(e-"‘r) 


where 


n j  =  min{<7,  p'  / 2} 


Hence  (6.32)  implies  that  the  asymptotic  relation 


7(r)  =  0(e“'^) 


holds,  which  concludes  the  proof.  | 


as  t  — ►  oo 


Proposition  6.1:  Under  the  assumptions  (I)  -  (TV),  the  system  (6.1T)r  has  a  unique  solution 
(«r(r),  /?r(r),  7r(r)),  for  r  >  0  and  r  =  Moreover,  there  exist  coefficients  ,  r  — 

1, . . . ,  N  of  the  expansion  (4.1a)  of  £*(e)  and  positive  numbers  ay  such  that 

°r(T)  =  0(e-ff'r)  'j 

/?r(T)  =  0(c"^r)  1  r  =  1 . JV  (6.35) 

7r(r)  =  0(c~"'r)  } 


as  t  — >  oo. 

Proof:  The  first  part  of  this  proposition  is  a  direct  consequence  of  Lemma  6.5  since,  by 
Assumption  (IV),  -Di/a^Mr))  is  nonsingular  for  all  r  >  0  and  by  the  derivation  of  the  system 
(6.17),  P0lr(T),  Pp,r(T)  and  P7iI.(t-)  are  well  defined  for  r  >  0  provided  only  the  previous  terms 
a*(r),  /?;(t ),  7,(r),  t  =  1, . . . ,  r  —  1,  are  available. 

Hence  it  remains  to  show  that  with  properly  chosen  coefficients  £*,  r  =  1,...,JV,  the 
asymptotic  relations  (6.35)  hold  for  all  0  <  r  <  N.  This  can  be  shown  inductively  by  applying 
the  methods  used  to  prove  the  second  part  of  Lemma  6.5.  Indeed,  Assumption  (IV)  ensures 
that  the  estimate  (6.35)  is  valid  for  r  =  0.  Thus  assume  that  there  exist  positive  numbers 
such  that  the  asymptotic  relations  (6.35)  hold  for  all  1  <  r  <  k  .  Note  that  the 
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exists.  Therefore,  by  integrating  tin:  first  equation  of  (6.17)*+j 


we  see  that 


<*Jt+i(r)  =  f/t+i  -  &+1  +  j  Pak+i(x)d*  ~  Pak+\(s)ds. 
Accordingly,  we  choose 

€t+i=6+l+  /  Pa,k+\(s)ds, 

Jo 


(6.38) 


whence 


a*+i(r)  =  -  J™  Pa,k+i(s)ds  =  -  J™  0(c~<^)ds  =  0(€~y^r) 


as  t 


oc 


and  the  asymptotic  relation  (6.35)  holds  for  ojt-n(r). 

Since  /?«(r),  70(r)  — +  0  as  r  — *  oo,  it  follows  that 

hrn^  {D^/affiofr))  -  DJ2(U0(r))(DJ3(fl0{T)))-'1 Dyh(il0(r))}  =  B( 0) 

where  by  Assumption  (III)  all  eigenvalues  of  £(0)  remain  strictly  in  the  left  half  plane.  Hence, 
by  applying  Lemma  6.5  to  (6.17)jt+i,  with  £J+!  chosen  by  (6.37),  we  obtain  the  existence  of 
ffJt+i  >  0  (  <  <7^,  )  such  that  the  other  two  components  /?*+1(r)  and  7*+i(t)  satisfy  the 
asymptotic  relations 

/?*+1(r)  =  0(e-^)l 

>  as  r  — »  oo. 

7fc+,(T)  =  0(e-'‘+*r)  J 

Thus  altogether  we  proved  that  the  existence  of  and  of  <7*+]  >  0  such  that  the  asymptotic 
relation  (6.35)  holds  for  r  =  k  +  1  and,  therefore,  (6.35)  is  valid  for  all  0  <  r  <  N.  | 


6.3.  The  procedure  for  generating  outer  and  inner  solutions.  The  proofs  of  Lemma 
6.5  and  Proposition  6.1  describe  the  procedure  for  the  generation  of  the  sequences  £*  and 

(Xr(t),V’r(0.Zr(0).  r  =  l . N  (6.39a) 

(ar(r),fl.(T),7r(r)),  r  =  l,...,Ar  (6.396) 

which  solve  the  systems  (4.2)r  and  (6.17)r,  respectively,  and  satisfy  the  asymptotic  relations 
(6.35)  provided  the  first  terms  (Xo(f),  ^o(0>  ^o(O)  and  (0,  /?o(r),  7o(T))  are  available  and 
satisfy  Assumptions  (I)  -  (IV).  Suppose  (Xr(t) ,  Yr(t),  Zr(t))  and  (or(r),  /?r(r),  7 At))  are 
available  for  all  r  <  k  and  satisfy  the  requirements.  Note  that  the  polynomial  P„,t+i(r)  in 
oi , . . .  ,7*  depends  only  on  the  values  of  (Xr(t),  Vr(f),  £»■(*))  and  their  derivatives  at  t  =  0  for 
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r  =  0, 1, . . .  ,k,  and  has  the  rate  Pa,k+i(T)  =  0(e-<7*+,r)  as  r  -+  oo.  Thus  Pa,fc+i(r)  is  known, 
and  furthermore  „ 


exists.  Then  with  specified  by  (6.38),  with  which  we  obtain  the  solution  a*+i(r)  of  the 
first  equation  of  (6.17)jt+i-  Prom  the  proof  of  Proposition  6.1,  it  follows  that  Ofc+1(r)  satisfies 
the  asymptotic  relation  (6.35).  With  the  initial  condition  Xfc+i(0)  =  £J+1  we  can  find  the 
unique  solution  (Xk+i(t),  Yit+i(i),  Z*+i(t))  of  the  system  (4.2)*+i.  Only  at  this  moment, 
Pffk+i(T)  and  P7fc+1(r)  are  fully  determined  because  they  depend  not  only  on  the  previous 
solutions  but  also  on  the  current  one,  (A*+i  ,lfc+i,  Zjt+i)-  By  inserting  ajt+i  =  ajt+i(T)  into 
the  last  two  equations  of  (6.17)*+i,  together  with  the  initial  conditions 

Pk+ l(O)  =  T)k+1  -  Pfc+i(0),  7*+i(°)  =  Cjh-1  -  Zk+ l(O) 

we  obtain  the  unique  solution  (/3fc+i(r),7t+i(r)).  In  other  words,  our  procedure  has  the 
following  form. 

Procedure  (A):  Let  (X0(t),Y0(t),  Z0(t))  and  (/9o(t)»7o(t))  be  the  solutions  of  the  re¬ 
duced  problem  (3.1)  and  of  the  problem  (6.17) o,  respectively,  which  satisfy  the  Assumptions 

(I)  -  (IV).  Forr  =  1 


(1)  solve  the  ODE 


to  obtain  or(r)  for 


*r(0)  =  fr-C 


=  tr  +  J  Pc,r(-S)ds; 


(2)  with  the  initial  condition  Xr(0)  —  £*  solve  (4.2)r  to  obtain  the  unique  solution  (Xr(t), 

Yr(t),  Zr(t)); 

(3)  insert  ar(r )  into  the  last  two  equations  of  (6. 17)r,  to  obtain  the  unique  solution  (/5r(r), 
7r(r)). 

Note  that  this  procedure  is  indeed  independent  of  the  concepts  of  outer  and  inner  solutions. 
But  the  sequences  £*  and  (6.39a, b)  generated  by  this  procedure  can  be  used  to  construct  outer 
and  inner  solutions.  Next  chapter  shows  that  inner  solutions  exist  and  can  be  expanded  with 
(ar(r),/?r(r),7r(r))  as  coefficients. 

7.  Existence  of  inner  solutions 

7.1.  Theorems  of  existences.  The  following  theorem  shows  that  the  series 
N  N  N 

J2ar(T)er,  X>(rK,  Y,  7r(ry 
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generated  by  Procedure  (A)  can  be  used  to  approximate  an  inner  solution  (a(r(e),  fi(r,  f), 
7(r,e)). 

Theorem  7.1:  Under  Assumptions  (I)  -  (IV)  suppose  that  the  sequences  £*,  („Yr(t),  Yr(t), 
Zr(t))  and  (ar(r),  fir(r),  7r(r)),  r  =  0, 1, . . . ,  N,  are  generated  by  Procedure  ( A ).  Let  (.Y(t,  e), 
Y(t,e),  Z(t.,e))  €  T>  {or  t  €  (0,  T]  be  the  outer  solution  guaranteed  to  exist  by  Theorem  5.1 
which  satisfies  the  initial  condition 

TOO=no~fo  +  £c.v, 

i=i 

and  the  asymptotic  relations  (5.1 ).  Then  there  exists  a  unique  solution  (ft(r,€),  0(t,  e),  7 (r,e)) 
for  the  inner  problem  (3.7)  on  the  interval  0  <  r  <  T/e,  for  all  0  <  e  <  fj ,  where  e,  (<  e0)  is 
a  sufficiently  small  number  such  that 

A"(t,f)  +  o(t/e,c)  €  Vx, 

Y(t,<)  +  fl(t/t,e)€  Vy, 

Z{U*)  +  7(t/e,e)  €  T>z, 

for  all  t  €  [0,  T]  and  0  <  e  <  ei ,  and 

N 

«(r-«)  -  53ar(T)fr  =  0(fw+1), 

r=0 

N 

/?(r,e)-£/?r(r)er  =  0(e"+1),  (7.1) 

r=0 

N 

l(r,f)-J2lr(ry  =  0(cN+'), 

r= 0 

uniformly  for  0  <  r  <  Tjc  as  e  — »  0. 

The  proof  of  this  theorem  will  be  given  after  next  theorem.  From  Theorem  5.1  and 
Theorem  7.1  follows  the  main  Theorem  of  this  paper. 

Theorem  7.2:  Let  Assumptions  (I)  -  (IV)  hold.  Then  there  exists  an  f  >  >  0.  0  <  f  j  <  e,, 
where  is  defined  in  Theorem  7.1,  such  that  for  all  0  <  e  <  f2  the  singularly  perturbed  DAE 
(l.la,b)  has  a  unique  solution  x  =  :r(t,  c ),  y  =  y(t,  e),  2  =  z(t,  c)  in  V  on  the  interval  0  <  t  <  T 
.  Moreover,  for  any  natural  number  N,  there  exist  an  outer  solution  (X(t,€),Y(t,e),Z(t,e)) 
€  V,  and  an  inner  solution  (a(r,  e),/?(r,  e),7(r,  c))  such  that 

x(t,  e)  =  X(t,e)  +  o  (t/e,  e),  y(t,  f)  =  Y(t,  c)  +  0{t/c,  f),  z(t,  f)  =  Z(t,  e)  +  7(t  A  ■  <0 
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for  0  <  t  <  T  and  0  <  e  <  cq.  Moreover,  the  outer  solution  satisfies 

N 

X(t,e)~Y,Xr(T)er  =  0(eN+1) 

r=0 

N 

Y(‘<  e)~J2  ¥r(TV  =  0^N+'  ) 

r=0 

Z((,£)-^Zr(r)e'  =  0(e"+>) 
r=0 

uniformly  for  0  <  t  <  T,  where  (Xr(t),Yr(t),  Zr(t))  are  determined  by  the  system  (4.2)r, 
r  =  0, 1, . . . ,  N,  and  the  inner  solution  satisfies 

N 

«((/«,€)  -'E«r(t/'W 

r=0 

N 

P(tJ*,*)-'ZMt/‘V  =  0(eN+') 

r=0 

N 

7(</«»  0  -  7r(</e)er  =  0(eN+1) 

r=0 

uniformly  for  0  <  t  <  T,  where  (afr(r))/?r(r),7r(r))  are  determined  by  (6.17)r,  and  satisfy 
° r(r)  =  0(e~efrT),  fir(r)  =  0(e~<r'r),  7 r(r)  =  0(e-ovr),  ,  r  =  0, 1, . . . , N,  as  r  ->  00 

for  some  or  >  0. 


7.2.  Proof  of  Theorem  7.1.  In  analogy  to  the  proof  of  Theorem  5.1,  the  first  step  here 
introduces  a  change  of  variables  into  the  inner  system  (3.7)  in  order  to  change  it  to  a  technically 
simpler  form.  Let  u  €  Rm,v  €  Rn,w  G  Rk  be  defined  by 

N 

a  =  u  -f  ^  or(r)er, 

r=0 

N 

0  =  v  +  Y^MtV +  A1{t)u,  (7.2) 

r=0 

N 

7  =  w  +  ^2  Jr(r)er  +  A2{t)u  +  Bi(r)v, 

r=0 

where  Ai(r ),  A2(t),  and  B\(t)  will  be  determined  such  that  after  the  change  of  variables  (7.2) 
the  system  (3.7)  will  simplify.  At  first  we  have  to  determine  again  a  proper  domain  for  the  new 


33 


variables  such  that  the  transformation  (7.2)  makes  sense  for  the  inner  problem  (3.7). 

For  this  recall  that  the  inner  system  (3.7)  was  obtained  by  introducing  the  change  of  variables 
(3.6a/b)  into  (1.1)  where  ( X{t,e),Y{t,* ),  Z(t,e))  is  an  outer  solution.  Thus  the  change  of 
variables  (7.2)  in  the  inner  system  (3.7)  is  equivalent  with  the  following  change  of  variables  in 
the  original  DAE  (1.1): 

N 

3-  =  .Y(t,e)  +  u  +  ^or(</e)€r, 
r— fl 
N 

y  =  Y{t,  e)  +  v  +  Y,  A-(*/0«r  +  M  (■ t/()u ,  (7-3) 

r=0 

N 

z  =  Z(t,e)  +  w  +  5^7r(</e)«r  +  A2(t/e)u  +  Bi(ife)v. 

r=() 

It  will  be  shown  that  there  is  a  small  neighborhood  Os  of  (0, 0, 0)  such  that  if  (u,  v,  w )  remains 
in  Os,  then  {x,y,z)  €  V  where  (x,y,z)  is  defined  by  (7.3).  For  this,  we  show  first  that  under 
Assumptions  (I)  -  (IV)  there  exists  a  sufficiently  small  0  <  f3  (  <  *2  )  RUCh  ^,at  Oo(0  +  A>(*  A)» 
Z0(t)  +  7 0(t/e))  eVyxVt  for  all  t  €  [0,T],  0  <  t  <  e3-  Indeed,  since,  by  Assumption  (IV), 
(Fo(0)  +  0o(t),  Zo(0)  +  7o(0)  €  Vy  x  2>,  for  all  r  >  0  and  V„  Vz  are  open  sets,  there  exists 
t*  G  (0,T)  such  that 

(yo(0  +  /?o(t/e)i^o(0  +  To(*A))  €  vy  x  for  all  t  €  [0,t*],0<  e  <  e3-  (7-4) 

On  the  other  hand,  since  (V0(f),  ^o(t))  £VyxVz  for  all  t  €  [0,T]  and  [0,T]  is  a  closed  set, 
there  exists  r  >  0,  independent  of  f,  such  that 

B(Y0(t),r)  x  B(Z0(t),r)  C  Vy  x  V„  V<  6  [0,T],  (7.5) 

From  limr_00(^o(7')i  7o(r))  =  (0,0),  we  find  that  there  exists  0  <  e'  (<  e3)  such  that 

\Po(t/e)\  <  r,  |7o(«/f)|<r,  Vt  €  [t*,r),0  <  e  <  4  (7-6) 

Hence  (7.5)  and  (7.6)  show  that 

(Xo(t)  +  Po{t/e),Z o(t)  +  7o(*A))  e  Vy  xV>'  for  a11  1  €  tr’Tl’0  <  e  (7‘7) 


which,  together  with  (7.4),  means  that 

(Y0(t)  +  p0(t/e),Z0{t)  +  7o (</0)  €  x  Vz,  for  all  *  6  (0,T],0  <  e  <  4- 


(7.8) 


By  Theorem  5.1,  the  estimates 

X(t,e)  =  X0(t)  +  O(t)\ 

y(t,e)  =  y0(t)  +  O(e)  1  as  e  — *  0,  (7.9) 

Z(t,e)  =  Z0(t)  +  O(e )  J 

hold  uniformly  in  t  G  [0,T].  Therefore,  (7.8),  together  with  (7.9),  implies  that  there  exists 
>  0  (<  £3)  such  that 


+  /30(t/e),Z(t,e)  +  y0{t/e))  e  Vx  xVy  xVz,  Vt  e  [0,t*],  0  <  e  <  e4.  (7.10) 

Note  that  (c*i(f/e), &(//«),  7i(i/e)),  i  =  1  are  uniformly  bounded  for  all  t  €  [0,T] 

and  0  <  e  <  e4.  Moreover,  suppose  that  A^(t),A-j(t)  and  Bx(r)  are  uniformly  bounded  for 
0  <  r  <  T/e.  Then  there  exist  65  >  0  (  <  e4  )  and  5  >  0  such  that  for  any  (u,v,iv)  6 
5(0,5)  x  5(0,5)  x  5(0,5),  the  point  ( x,y,z )  defined  by  the  transformation  (7.3)  is  contained 
in  V  for  all  t  €  [0,  T]  and  0  <  e  <  65.  We  write  Pu  =  -#(0,5)  C  Rm,  Vv  =  5(0,5)  C  5", 
T>w  =  5(0, 5)  C  21* ,  and  V  =  Vu  x  Z>„  x  T>w. 

The  change  of  variable  (7.2)  is  now  introduced  in  (3.7).  For  abbreviation,  the  notations 
fti(r),  r,(r)  defined  in  (6.10)  and  (6.11)  will  be  used  again.  Furthermore,  set 


A  = 


Then 


/  / 

0 

0 

0\ 

/«\ 

Ai{r) 

I 

0 

0 1 

\  Mr) 

Bi(r) 

I 

or 

M 

V  0 

0 

0 

1/ 

Vo/ 

n=  y  =n0(r)  +  ^f2<(r)ei  +  A(t)U  +  Rn(r,e) 


and 


r  = 


^(fr’c)\  N 


(7.11) 


(7.12) 


where  5„(r,  e)  =  0(eiV+1)  is  independent  of  (u,v,w).  Thus,  by  expanding  /(ft)  and  /(r)  at 
fto(r)  and  r0,  respectively,  we  obtain 


/(ft)-/(r) 

*-l  3' 

>  Pili To) 


J=1 
N 

J=1 


^  ft*(r)e*  +  ,4(r)2/  j  +  5(r,u,  v,  u>,  e) 

*-l 

Er*(r)<‘)  -»( r,e). 


(7.13) 
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As  a  remainder  term,  i?(r,  w,  i?,Ti>,f)  satisfies 

(i)  i?(</«,0,0,  0,  «)  =  0(tN),  as  e  — ►  0  (7.14a); 

(ii) 

DuR(t/(,u,v,w,e)  -  0(e  +  |tz|  +  |f|  +  |te|)  > 

DvR(t/f, «,  V,  Wj  e)  =  0(t  +  |u|  4-  |e|  +  |u>|)  >  as  e,  |u|,  |i>|,  |u>|  — »  0  (7.146) 

DwR(t/e,u,v,  tv,  e)  =  0(e  +  |u|  +  |u|  +  |u>|)  J 

uniformly  for  0  <  t  <  T. 

Moreover  the  remainder  term  H(t,(.)  =  0(eN+1)  satisfies  (7.14a,b)  as  well.  Obviously,  the 
conditions  (7.14a,b)  are  exactly  the  same  as  Condition  (N).  It  follows  from  (7.13)  that 


rn)  -  /( r) 

N 

=  /(Ho(r))  -  /(r0)  +  £crer  +  Df(Q0(T))A(r)U  +  G(r,u,r,te,  t) 


where  cr  is  as  defined  in  (6.15)  and 

G(t,  u,  t;,  w,  e)  =  J2(r,  u,  v,  te,  e)  +  0(eN+1 )  -  //(r,  e) 


}= 2 


DjW0(T)) 


/  N 


\i=i  4 


Under  our  assumption  that  A(t)  is  uniformly  bounded  on  0  <  r  <  T/(,  we  see  that 


N 


E 


Djmo(r)) 

j'- 


j>*(r)«*  +  A(t)U 


) 


satisfies  (7.14a, b).  Since  the  Condition  (N)  also  holds  for  R(r,u,v,  «*,  e),  0(t^+1),  H(t,  e), 
this  implies  that  G(t,  u,  v,u>,e)  satisfies  Condition  (N)  as  well. 

Substituting  (7.2)  into  the  first  equation  of  (3.7a)  we  obtain 


Tr+t  ^ -•(/,(«)  ~MD) 

(7.16) 

=  e[/i(fto(r))  -  /i(r0)  +  £rrer  +  D/,(O0(r))^]  +  G,(r,  mr,  ,e,  f). 

r=l 

Since  «,(r)  satisfies  the  first  equation  of  (G.17)r,  it  follows  from  (7.16)  that 

~  =  tDJ\{Sl0{T))A{T)U  +  C7i(r,u,t;,u?,e)  (7.17) 
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where  G\  satisfies  (7.14a, b).  Similarly,  substituting  (7.2)  into  the  second  equation  in  (3.7a), 
and  utilizing  the  fact  that  /3r(r )  is  the  solution  of  (6.17)r,  r  =  1, . . . ,  N,  we  see  that 

dv  . »  >  \  ,  , 

-+Mt)u  +  Mt)-  (718) 

=  D/2(S2o(t)).4(t)W  +  G(2)(e,  u,  u,  w,  r) 

where  G(2)is  defined  by  (7.15)  with  /  replaced  by  /2.  Obviously,  the  function  on  the  right  side 
of  (7.18)  satisfies  (7.14a, b).  Thus 

G2(e,u,v,w,T ) 

=  -eA1(r)Df1(il0(T))A(T)U  -  Ai(t)Gi(t,u,v,w,  e)  +  G(2)(e,u,u,u;,r) 
satisfies  (7.14a, b)  as  well,  and  it  follows  from  (7.18)  that 

=  Df2(£lo(r))A(T)H  -  Ai(r)u  +  G2(e,u,v,w,T).  (7.20) 

dr 

Finally,  substituting  (7.2)  into  the  third  equation  of  (3.7a),  we  find  that 

DMS10{t))A(t)U  +  G3(r,u,  v,  w,e)  =  0  (7.21) 

where  G 3  is  defined  by  (7.15)  with  f  replaced  by  f3. 

Note  that  from  the  derivations  of  the  equations  (7.17),  (7.20)  and  (7.21),  we  see  that  the 
functions  Gj,  i  —  1, 2,3,  are  well  defined  on  V  x  (0,  e5]-  Now  we  show  that  there  exist  matrices 
Ai(r),  A2(t)  and  Bj(r)  which  are  uniformly  bounded  on  0  <  r  <  T/e  and  for  which  the 
equations  (7.17),  (7.20)  and  (7.21)  simplify  considerably.  For  this,  we  write 

Df2{Sh{r))A{T)U  -  A\{t)u 

=  [Dxf2(n0(r))  +  DyMSloirVA^T)  +  DJ2(S10(t)A2(t)  -  Aj(r)]u 
+  [Dyf2{n0(T))  +  D,/2(fto(r))B,(r)]v  +  DJ2(fl0(r))w 


Df3(Cl0(r))A(r)U 

=  [B,/s(fto(r))  +  By/3(fto(r))A,(r)  +  D*/3(fto(r)A2(r)]u 
+  [J3,/j(fio(r))  +  DJ3(^o{t))Bi(t)]v  +  DJ3(^o(t))w, 
and  use  the  equations 

B,/2(n0(r))  +  B,/2(fto(r))i4,(r)  +  B,/3(i2o(r))A2(r)  -  A\(r)  =  0, 
B«/s(fto(r))  +  Dy/3(ft0(r))A,(r)  +  D,/3(fio(r))A2(r)  =  0, 
Dsf3(Sl0(r))  +  B4/3(f20(r))B1(r)  =  0, 
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to  determine  -Ai(t),  .^(t)  and  Bj(r).  In  order  to  show  that  the  solutions  ^(t),  A2(t)  and 
Bj(t)  of  (7.23)  are  uniformly  bounded  for  0  <  t  <  T/e,  note  first  that  the  third  equation  of 
(7.23)  gives 

B,(r)  =  -(DJ3(Sl  o(r)))-1  Dyf3(Q0(r)) 

which  is  obviously  bounded  for  0  <  r  <  T/t.  Indeed,  there  is  a  constant  C,  independent  of  r 
and  e,  such  that  ||Bi(r)||  <  C  for  all  r  >  0.  The  second  equation  of  (7.23)  now  provides  that 

A2(t)  =  -  (D,f 3(Q0(T)))~i  (DzMQ0(t))  +  DyMClo(r))A,(T)). 

By  substituting  this  into  the  first  equation  of  (7.23)  we  obtain  that 

^ii,(r)  =  Af(r)A,(r)  +  C(r) 

where 

M(t)  =  Dyf2(il0(r))  -  DtMQ0(r))(DJ3mr)))  By/3(f20(r)), 

C(r)  =  Dzf2(il0(r ))  -  D,f2(Q0(r))  (Dzf3(Q0(r)))  Dtfs(Si0(r)), 

Because  of 

lim^Af(r)  =  B(  0) 

and  5?(A(B(())))  <  0  for  all  eigenvalues  A(B(0))  of  B( 0),  tliere  exists  a  r0  >  0  such  that  all 
eigenvalues  of  M(r)  remain  strictly  in  left  half  plane  for  r  >  t0.  Thus  if  -4j(r)  is  chosen  as  the 
solution  of  (7.25)  that,  satisfies  the  initial  condition  >h(0)  =  7,  then  Lemma  5.1  ensures  that 
A](t)  is  uniformly  bounded  in  r,  t  >  0.  With  this  Aj(r)  it  follows  directly  from  (7.24)  that 
A2(t)  is  uniformly  bounded  in  t  as  well. 

Note  that 

N 

u(0,  e)  =  o(0,  e)  -  ^  ar(0)er 

r=0 

N 

=  no -£«.>- cv 

r=0 

=  6^)  =  ()(<"+'),  (7.27a) 

N 

i;(0,  f  )  =  /?( 0,  f)  -  £  Prioy  -  i4i(0)«(0, 0 

r=0 

=  02(t)  =  O(eN+'),  (7.276) 

N 

™(0,  f )  =  7(0,  e)  -  ^  7r(0)f r  -  i42(0)«(0,  f)  -  B]  (0)e(0,  f ) 

r=0 

=  9*{*)  =  0{tN+1).  (7.27c) 


(7.24) 

(7.25) 

(7.26) 
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Therefore,  with  these  specially  chosen  A](t),  A2(r)  and  Bi(t),  it  follows  from  (7.17),  (7.20), 
(7.21)  and  (7.27a, b,c)  that  (3.7)  becomes 

^  =  e(Ci(r)u  +  Ti(r)w  +  J5i(r)u>)  +  eGi(e,  u,  v,w,t) 
dr 

—  =  M(t)v  +  E2(t)w  +  G2(e,  u,  v ,  w,  t) 
dr 

0  =  E3(t)w  +  G3(e,  u,v,w,t)  (7.28) 

u(0)  e)  =  0i(e)  =  0(ew+1), 
v{0  ,e)  =  62(e)  =  0(eN+'), 
u>(O,e)  =  03(e)  =  O(eN+1). 

This  is  equivalent  to  the  following  integral  equation 

u(r, e)  =  $(r,0, e)6*i(e)  +  f  $(r,s,e)e(Ti(s)v(s, e)  +  £i(s)uj(s,e) 

Jo 

+  Gi(e,  u(s,  e),  u(s,  e),  u>(s,  e),  s))ds, 

v(t,  e)  =  9(r)ff2(e)  +  ¥(r)  f  (s)(E2(s)w(s,  e)  (7'29) 

Jo 

+  G2(e,  u(s,  e),  v(s,  e),  tu(s,  e),  s))ds , 
w(t,  e)  =  -257Hr)G3(e,u(r,e),v(r,e),to(T,e),r), 

where 

Ci(r)  =  DMCloir))  +  -Dy/iWO^r)  +  D,/i(flo(r))j42(r), 

Ti(r)  =  DM<l0{r))  + 

Ei(T)  =  DJ1(Cl0(r)),  E2(r)  =  DJ2(nQ(T)),  Es(r)  =  D2f3(Q0(r)), 
and  $(r,  $,e)  satisfies 


d$(r,  s,e) 


=  eCi(r)$(T,  s,e) 


dr  v  ’  ’  7  (7.30) 

$(r,s,e)|r=,  =  I 

while  ^(r)  is  a  solution  of  following  system 

f  =  (7.31) 

¥(0)  =  1. 

It  turns  out  that  $(t,s,  e),  ^(r)  are  uniformly  bounded  for  0  <  s  <  r  <  T/e  and  0  <  e  <  Co¬ 
in  deed,  because  Ci(r)  is  uniformly  bounded  for  all  t  >  0,  there  exists  a  constant  C  >  0  such 
that  |Ci(t)|  <  C  for  all  r  >  0.  Thus  from 


$(r,s,e)  =  J+ejf  Ci(A)$(A,s,e)dA 
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it  follows  that 

||*(r,  s,  e)||  <  1  +  eC  j*  ||$(A,  .s,  c)|| d\ 
which  by  Gronwall’s  inequality  implies  that 

||*(r,a,e)||  <  e'c(r_,)  <  eCT 

for  all  0  <  a  <  t  <  T/e.  Hence  $(7-,  s,e)  is  uniformly  bounded  for  0  <  s  <  r  <  T/e. 
For  the  matrix  M(r),  there  exists  a  positive  number  fi  such  that 


R(A(M(r)))  <  -//,  r>r0 

holds  for  all  eigenvalues  of  M(t).  Hence  by  Lemma  5.1  (  with  e  =  1  )  it  follows  that 
||'t(r)(4'(.s))-I||  <  Ke-^-W 

holds  for  0  <  s  <  r,  where  K  is  a  constant  which  is  independent  of  s,r. 

Now  we  show  that  (7.29)  has  a  unique  solution  (n,v,tv)  for  which 

ti(r1«)  =  0(ew+l),  v{r,f)  =  0(tN+1),  u>(t,  e)  =  Ofe"41 ), 

uniformly  for  0  <  r  <  T /e. 

With  the  change  of  variables 

u(M)  =  u(t/e,  e), 

=  v(t/e,e), 
w(t,  e)  =  w{t/( ,  e), 
t  ~  t/f,  a  =  s'/e, 

the  constrained  integral  equation  (7.29)  becomes 


■s: 


«(M)  =  ^(^0,  e)^i(c)  +  /  $(t,  Ae)(Li(s7<0£’(.s',e)  +  Ei(s' /()w{s\f) 


where 


+  G]  (e,  u(s\  e),  t>(s\  e),  te(s\  e),  s'))ds' 

*(<,  0  =  *(*,  t)62(€)  +  f  '(s’'€)(E2(s'le)w(s’,  e) 

Jo  € 

+  G2{f,  «(*\  f),  (),  w(s',  ( ),  s'))ds' 

«(t,e)  =  -Ef1(t/e)G3(f,u(t,e),v(t,e),w(t,e),t) 


$(t,s',c)  =  3>(</e,  s/e,  e) 

fc^e)  =  *(t/f) 

Gi(f,u,v,w,t)  =  Gi(f,u,v,w,t/t),  for  i  =  1,2,3. 


(7.32) 


(7.33) 


(7.34) 


(7.35) 
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(7.36) 


Because  of  (7.33)  the  kernel  in  the  second  equation  in  (7.35)  satisfies 

||®(<,e)4rV,e)||  < 

for  all  0  <  s'  <  t  <  T.  Moreover  4(f ,  s',  e),  L^s'/e),  E^s'/e),  i  =  1, 2, 3,  are  uniformly  bounded 
in  0  <  s'  <  t  <  T.  Since  Gj, »  =  1,2,3,  satisfy  (7.14a, b),  we  see  that  Gi(t,u,v,w,t),  i  =  1,2,3, 
satisfy  the  Condition  (N).  Now,  with  the  same  notation  as  in  (5.15)  we  rewrite  (7.25)  in 
the  form  of  constrained  systems  of  integral  equations  as  the  one  in  (Ya2).  Then  Theorem  1,2 
in  [Ya2]  guarantee  that  the  system  (7.35)  has  a  unique  solution  (u(<,  e),u(t,e),u>(<,  e))  on  the 
interval  0  <  t  <  T  for  which 

u(M)  =  0(ew+1)  ] 

6(<,e)  =  0(eN+1)  1  as  e  — >  0, 

w(t,e)  =  0(eN+1)  J 

uniformly  on  0  <  t  <  T.  Hence  the  system  (7.29)  has  a  unique  solution  (u(r,e),t;(r,e),u>(r,  e)) 
on  the  interval  0  <  r  <  T/e  which  satisfies  the  estimate  (7.34).  This  completes  the  proof  of 
Theorem  7.1.  g 
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